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Abstract 



We compute the QCD next-to-leading order matching conditions of the (chromo)- 
magnetic operators relevant for B — > X s ^y in super symmetric models with mini- 
mal flavour violation. The calculation is performed under the assumption that the 
charginos and one stop are lighter than all other squarks and the gluino. In the 
parameter region where a light charged Higgs boson is consistent with measure- 
ments of BR(5 — > X a j), we find sizeable corrections to the Wilson coefficients. 
As a consequence, there is a significant reduction of the stop-chargino mass region 
where the supersymmetric contribution has a large destructive interference with the 
charged-Higgs boson contribution. 



*On leave of absence from INFN, Sez. di Padova, Italy. 



1 Introduction 



The inclusive decay rate for B — > X s7 has been first measured by CLEO with 
the result BR(B -> X s7 ) = (2.32 ± 0.57 stat ± 0.35 syst ) x KT 4 [|. Recently, a 
preliminary new result based on about 30% more data has been presented by the 
collaboration, BR(B -> X s7 ) = (2.50 ± 0.47 stat ± 0.39 syst ) x 10~ 4 @. The same 
process has also been measured by ALEPH at LEP, with the result BR(i? — > X s7 ) = 
(3.11 ±0.80 

stat 

±0.72 syst ) x 10~ 4 1. 
There has been significant theoretical effort in refining the prediction of BR(B — > 
Xg'y) in the Standard Model (SM). Calculations are now available for the next-to- 
leading order (NLO) corrections to the anomalous dimensions 0, the matrix el- 
ements ||, and the matching conditions of the Wilson coefficients |6|, [7|, for the 



leading non-perturbative effects ||, and for the QED corrections || pTOH . In par- 
ticular the QCD NLO corrections are important since they reduce the large scale 
dependence of the leading result jlTf , which amount to a 30% uncertainty in the 
leading order (LO) theoretical prediction. Recently, Kagan and Neubert [|Kj have 
argued that, at the NLO, the contributions to BR,(B — ► X s7 ) from the different 
Wilson coefficients exhibit a larger scale dependence than the total result, signal- 
ing that the theoretical uncertainty may have been underestimated in previous lit- 
erature J|, [12], 0| . Nevertheless, even considering this effect, they find that the 
theoretical error can be evaluated to be about 10%. Combining the different the- 
oretical studies, the recent complete analysis in ref. JED] gives the SM prediction 
BR(B -> X s7 ) = (3.29 ± 0.33) x 10" 4 x BR(B -> X c eu)/10.5%. 

While the theoretical prediction is made for the total rate of BR(B — ► X s j), 
the experimental data refer only to events with photon energies between 2.2 and 2.7 
GeV. The extrapolation from the data to the total rate introduces further theoretical 
uncertainties from the calculation of the photon energy spectrum [13]. This aspect 
has been recently emphasized by Kagan and Neubert [l(| who have recomputed the 
spectrum and extrapolated from the CLEO high-energy photon data the total rate 
BR(B -> X s7 ) = (2.66 ± 0.56explo:« th) x !0" 4 - This is compatible with the SM 
within one standard deviation. 

To summarize, although the SM prediction is still higher than the CLEO result, 
the discrepancy seems no longer statistically significant. We should however wait 
for the forthcoming CLEO analysis and future studies at LEP to further clarify 
the situation. Meanwhile, we use the conservative CLEO upper limit of BR,(B —> 
X s ~f) < 4.2 x 10 -4 at 95% C.L. in our numerical analysis. 



1 



The inclusive decay B — > X s ^ is a particularly interesting probe of physics 
beyond the SM, since it is determined by a flavour-violating loop diagram. Because 
of this welcome sensitivity on new physics, it is important to determine the predicted 
rate with sufficient accuracy in a variety of models. In the case of two-Higgs doublet 
models, results for the next-to-leading order (NLO) QCD corrections have been 
presented in refs. J7J, |TJ], [TJ|. In these models, a new contribution arises from charged- 
Higgs exchange which always increases the SM prediction for BK(B — > -XVy). As a 
consequence, strong limits on the charged-Higgs mass can be derived. These limits 
are quite dependent on the treatment of theoretical errors and on the amount of 
discrepancy between theoretical and experimental results [[7], . 



Here we want to extend our previous analysis of QCD NLO effects in BK(B — > 
X s ^f) to the case of supersymmetry. Refining the theoretical calculation in supersym- 
metry [16, [l7j (for reviews see ref. [18]|) is important because BR(5 — > X s ^) provides 
a very stringent constraint on the parameter space of the model. Moreover in super- 
symmetry it is possible to evade the strong limits on the charged Higgs mass, since 



the chargino contribution can interfere destructively and reduce BK(B — > X s ^) [IT? 



Our result allows us to give a more reliable estimate of the model parameters nec- 
essary to achieve this cancellation. 

We have not tried to perform a complete analysis in the most general supersym- 
metric model, an effort, we feel, that is penalized by the fact that the final result 
will be very involved and will depend upon many unknown parameters. Instead, we 
try to focus on what we believe is the most interesting part of the parameter space 
for B — > Xg'-f. Specifically, we concentrate on the case in which the flavour violation 
is completely dictated by the Cabibbo-Kobayashi-Maskawa (CKM) angles (insuring 
predictability) and in which the charginos and a scalar partner of the top are light 
(insuring a sizeable new contribution). 

The paper is organized as follows. In sect. 2 we define our assumptions on 
the supersymmetric parameters and discuss the underlying hypotheses. Section 3 
contains the analytical results for the supersymmetric corrections to the relevant 
Wilson coefficients at the matching scale. A numerical analysis of BK(B — > X^) in 
supersymmetry at the NLO is presented in sect. 4. 
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2 Minimal Flavour Violation and Heavy Squark- 
Gluino Effective Theory 

In this section we describe our main theoretical assumptions on the supersymmetric 
model: i) the minimal flavour violation and ii) the heavy squark-gluino effective 
theory that describes chargino, stop, and the SM degrees of freedom with other 
supersymmetric particles integrated out. 

In the supersymmetric extension of the SM with general soft-breaking terms, 
there are a variety of new sources of flavour violation. The supersymmetry-breaking 
squark masses and trilinear terms lead to a mismatch in flavour space between 
quark and squark mass eigenstates. The flavour violation is then described by a 
large number of mixing angles, which cannot be determined theoretically. All we 
know is that present measurements and limits on various flavour- violating processes 
provide quite stringent bounds on these mixing angles (see e.g. ref. |19||). 

A more predictive case is what we will call minimal flavour violation, in which 
all flavour transitions occur only in the charged- current sector and are determined 
by the known CKM mixing angles. This is indeed the case in several theoretical 
schemes in which the communication of the original supersymmetry breaking to the 
observable particles occurs via flavour-independent interactions. In many of these 
schemes (especially when supersymmetry breaking occurs at low-energy scales) the 
departure from the minimal flavour violation hypothesis, caused by quantum effects, 
is rather small. Our assumption can then be justified in gauge-mediated models p0 | 



(for a review, see ref. [21|) and in certain classes of supergravity theories. At any 
rate, it plausibly corresponds to the unavoidable flavour violation, present in any 
supersymmetric model. 

Within this class of models, we focus on the special case of interest for B — > X s ~f, 
in which one stop is considerably lighter than the other squarks. In order to pre- 
serve the successful fit of the electroweak precision measurements, we assume that 
the light stop is predominantly right-handed P^fl . This assumption can also be the- 
oretically justified by the observation that the renormalization-group evolution of 
the squark mass parameters indeed pushes the right-handed stop mass to smaller 
values. Therefore, we concentrate on the case in which the supersymmetry-breaking 
parameters in the up-squark sector are flavour diagonal (but not necessarily univer- 
sal) in the basis in which the corresponding up-type quark mass matrix is diagonal. 
The flavour-violating stop interactions arise only from charged-current effects and 
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are completely determined by the CKM angles. The stop mass matrix is given by 

? / m~ L + m\ + (| - | sin 2 9 W ) cos 2/?M| m t (A - /i cot (3) \ 

m * ~~ V wit (A t - //cot/3) m 2 + m 2 + | sin 2 6 W cos 2/3Mf / 

Here m \ L -> an d At are supersymmetry-breaking parameters, fi is the Higgs 

mixing mass and tan /3 is the ratio between the two Higgs vacuum expectation values. 
We also assume that the charginos are lighter than gluinos and heavy squarks. These 
approximations considerably simplify the calculation and the final expressions, and 
they are appropriate to identify the leading supersymmetric contribution. 
To implement this scenario we assume the following mass hierarchy 

fig ~ Q{mg,mq,m^ > fi w ~ 0(M w ,m t ,m x ±,m i2 ) > m b > A QC d- (2) 

Here and in the following fi w is the mass scale of the charginos {x^), and of the 
lighter stop (t 2 ), to be identified with the ordinary electroweak scale. The charged 
Higgs (H^) mass is not constrained. The scale fig, characteristic of all other strongly- 
interacting supersymmetric particles (squarks and gluinos), is assumed to be larger, 
say of the order of the TeV. We compute the QCD NLO corrections to the Wilson 
coefficients keeping only the first order in an expansion in fi w j fi- g . 

As usual, the presence of different mass scales allows us to use a stack of effective 
theories, obtained by integrating out of the theory, at each matching scale, the heavy 
degrees of freedom. The effective theory just below the scale fi g is particularly 
simple. The only dimension-five operators involving supersymmetric particles are 
the chargino (chromo)-magnetic dipole moment and the operator 1 1 1\ t 2 , obtained 
by integrating out the gluino. However, an explicit calculation shows that these 
operators do not contribute to the NLO matching conditions of magnetic-dipole 
operators at the scale fi w . Therefore, to our end, the inclusion of the leading fi w j fi g 
corrections does not require any new operator other than renormalizable ones. 

The next step is the running of the intermediate effective Hamiltonian between 
fig and fi w . Although the relevant anomalous dimensions are known |23|], log re- 
summation is not really needed, given the smallness of the relevant log terms, 
a s(/ i w )/4vrlog(/i?//i 2 v ) ~ 0.01-0.05 for typical values of the supersymmetric masses. 

Non-resummed NLO effects of heavy supersymmetric particles, including I /fig 
corrections, are simply given by the one-loop Feynman diagrams, containing these 
particles in the loop, which contribute to renormalizable operators in the intermedi- 
ate effective theory. In other words, we need to compute only corrections involving 
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gluinos and heavy squarks to the coupling constants appearing in the vertices in- 
volving the chargino, x^M^i ^ ne ^ vector boson, W~qq', the charged physical, 
H~qq', and unphysical, <f)~qq', scalars. Heavy O(fig) particle effects renormalize 
the masses of 0(fi w ) particles, but these effects are reabsorbed through the defi- 
nition of renormalized masses. We use on-shell masses for squarks. Quark masses 
are also defined on-shell, except corrections not involving supersymmetric particles, 
which are subtracted in the MS scheme. This definition simplifies the insertion of 
supersymmetric contributions into existing RG-evolution formulae, where MS run- 
ning quark masses are commonly used. Because of our assumption that the light 
stop is mainly right-handed, the stop mixing angle % is small. For instance, tak- 
ing the supersymmetry-breaking parameter A t in eq. (|1|) of the order of fi g , yields 
Ql ~ 0(n w /n g ). In this framework, l/fi g corrections multiplied by sin^- factors 
should be regarded as of higher order. 

Results for the renormalized vertices in the intermediate effective theory are col- 
lected in appendix A. Notice that the Xa ^2 vertex renormalization is infinite. This 
is not surprising, since supersymmetry and the GIM cancellation are spoiled in the 
intermediate effective theory, where the heavy squarks and gluinos have been inte- 
grated out. This divergence cancels out in the matching of the magnetic operators 
at the scale fi w against a corresponding divergence generated by the insertion of the 
chargino vertex into two-loop diagrams. This cancellation actually provides a check 
for the results of the two-loop calculation. However, the presence of the divergence 
calls for a regularization. We choose the naive dimensional regularization (NDR), to 
be consistent with the calculation of the anomalous dimension matrix |3J . However, 
it is known that NDR breaks supersymmetry. In particular, the renormalizations 
in NDR of the gauge boson and gaugino interactions with matter, as well as the 
Higgs boson and higgsino interactions, are different and manifestly violate super- 
symmetry. Supersymmetric Ward identities are restored with appropriate shifts of 
the gauge and Yukawa couplings in the the x^M"2 vertex [p4j| , denoted as r\ Y and 
i] g in appendix A. These shifts correspond to the difference of using NDR versus 
dimensional reduction (DR) ||25|| , a regularization that preserves supersymmetry. 

The formulae in appendix A are given in terms of few functions of the ra- 
tios m~/m~. Notice however that, whenever m~ 2 /m~ appears, only terms up to 
0{mi 2 /rrig) should be retained in the corresponding functions to be consistent with 
the operator product expansion. Finally, we also keep the leading contribution in 
the sbottom mixing angle, since the left-right mixing term {A^ — /ztan/3) rrib becomes 
important for large values of tan f3. 
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3 Wilson Coefficients 



This section contains the result for the NLO supersymmetric contributions to the 
Wilson coefficients C** (pw)- Following the discussion of the previous section, the 
calculation of the NLO corrections to (nw) an d Cg (t^w) can be divided in two 
parts: i) the contribution of the heavy particles at the scale fj,g that will appear as 
renormalization of the coupling constants in the LO diagrams, ii) The contribution 
of the intermediate-scale particles that requires the computation of the two-loop 
gluonic corrections to the LO supersymmetric diagrams involving the charginos and 
the light stop. Concerning the latter two strategies are at hand. One can match 
matrix elements of operators belonging to a basis obtained enforcing the equations 
of motion, a procedure that however requires an asymptotic expansion of the rele- 
vant diagrams in the external momenta. Alternatively, one can use a larger off-shell 
operator basis and perform the matching on off-shell matrix elements. In this case, 
one can use the freedom of the off-shell status to choose a suitable kinematical 
configuration such that the various Feynman diagrams can be evaluated using or- 
dinary Taylor expansions in the external momenta. This second strategy, already 
successfully applied by us in ref. JjJ to the calculation of the QCD corrections to the 
matching conditions of the AB = 1 magnetic and chromo-magnetic operators in the 
SM and in two-Higgs doublet models, has also been employed in this calculation. 
Concerning the technical details we refer to ref. 0. 

In the supersymmetric model under consideration we can organize the Wilson 
coefficients of the operators entering the effective Hamiltonian in the following way: 

Ct"M = cl 0)eff ^ w )+5 H C^ ff ^ w ) + 5 s C^ ff ^ w ) 

+ fe (1)e// (M + s*cP*"M + *M 1W/ M (3) 

where Ci^(fj, w ) represents the SM contribution (k = 0,1), 8 H C^^{pL w ) the 
additional terms present in a two-Higgs doublet model, while w ) contains 

the contribution from supersymmetric particles. Explicit expressions for Cf )e// (/i 
and 5 H Cl h)ef, (ji w ) can be found in sections 4 and 5 of ref. [J7J 1 . Concerning the 
remaining contributions, the LO 5 s Cf (fiw) terms are given by |L6| 

8 s Cf )eff {yL w ) = i = l,...,6 



^J° 8 )e// (M = E 

3=1,2 



2M WT>2 F (1)/ \ | Uj2 M Wy p (3)/ N 



x In ref. 5 H c\ k)eff are indicated as 5C Weff 
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Ml 



r lj m 2 T 



2 



u. 



J2 



Mr, 



■tyCOSflfF&^J/y) 



V2cosf3 



rn 



Xj 



3 -2j ^X^7,8 12/2J, 



t 



12 



^2 M w 

\/2 cos (3 m 



t 2j sin 9tFj 3 i(y: 



2j) 



Xj 



(4) 



where the functions Fyg are defined in eqs. (29)-(30) of ref. and 



i(3), 



F 8 (3) (x) 



x 



5 — 7x x (3 x — 2) 

ro H : — in x 



We have defined 



£y = K U COS ^ 



m 

2 ' 



6(x-i; 

1 + X 

2(x-V 



tan 0{Yj] 



m~ 
tk . 

mi 



3 (x - 
x 



v. 



lnx. 



Y; 



V j2 cos 6 ] i m t (jj w ) 



(5) 

(6) 
(7) 



v^sin/? M M 

where m t (fi w ) is the top-quark running mass at the scale fi w . 

The stop eigenstates t\ = cos 9ft l + sin 9ft r and t 2 = — sin 9ft l + cos 6^ i R have 
mass eigenvalues and m f ~ 2 and we have taken all other squarks to be degenerate 
with mass m. The two matrices U and V diagonalize the chargino mass matrix 
according to (M is the weak gaugino mass) 



U 



M 

M w y/2cosp 



M w >/2sm(3* 



V 



-i 



and henceforth the scalar quark 



and are assumed to be real. Notice that in eq. 
masses are understood as on-shell. 

The last ingredients needed for a complete NLO calculation in our supersym- 
metric model are 5 s C\ l)eff ^ w ). For the current-current and penguins operators we 
have 



0* 



(_/ 4 \H-w j 





Ml 



1,2,3,5,6 



*2 3=1,2 



E t 2 2j E x (y: 



with 26 



E x (x) 



x (11 - 7x + 2x 2 



x 



X 



3(x- 1) 



2JJ 



lnx. 



(9) 
(10) 

(11) 
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The 0(a s ) corrections to the coefficient of the (chromo-) magnetic operator can be 
divided into 4 pieces 



°7,8 



+ ^<8 )e// (M 



(12) 



Here 5 x C^l e ^ represents the chargino contribution while (5 W , d^ 1 , 5^ 2 )C^ e ^ take 
into account the renormalization effects due to the O(fig) heavy particles in the W, 
physical and unphysical charged scalar couplings, respectively. We find 



5 x a 



(l)e// 



5 x cf e// (M 



3=1,2 



a M w 



Xj 



\/2 cos f3m Xj \ - in 
\Y 3 R s F?\y 2j )) - lY 3 ^ft 2j (R s + R b ) F?\y 2j ) 

*2 



6xcP eff (» w )[7^8,-~E x (y 2j ) 
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(13) 
(14) 



Eq. (|T4]) means that 5 X C$^ (j^ w ) can be obtained from the r.h.s. of eq. ( |T3"D by 
replacing the index 7 with the index 8 in the various functions and by substituting 
the term —4/9 E x (y 2 j) with — 1/6 E x (y 2 j). In eqs. (|i"3|)-([ll]) we have used 

Gf\x) = 4F 7 (1) (x)fAW+A«-l + 31n m ^ 



9 V m; , 

x (85 - 347a; + 526 x 2 ) 4x 2 (-8 + 13x + 6x 2 ) T ( 1 

+ ~ L + — — -Li 2 l-~ 

V x 



243 (1 - x) 



M{x-l) q 



9(x-l) 4 
lnx + 



4x(20- 126x+ 144x 2 + 39x 3 ) , 2 x 2 (21x - 10) , 2 , , 

- -^Wx (15) 



9(x-iy 



G^(x) 



+ 



,(3), 



x) A 



4:i 



A« + A® 



2) _16 (3 ■■rlr [A J i I 



9(x- 1)' 



4(3:r- 5) 4 (4 -30 x + 40 x 2 ) 16 (1 - 3x) x , 



9 (ar — 1)' 



9(x- l) c 



In a; — 



9 (x - 1)*' 



In x 



(16) 



-F^fx) 



AS + A« 



3 In 
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x (1210 -437s- 1427a; 2 ) x 2 (49 + 46 x + 9 x 2 

648 (x - l) 3 
x (85 - 603 x - 387 x 2 + 78 x 3 ) 



108 (x - 1)" 



lnx 



12(x- l) 4 
13 x 2 



Li 2 1 - - 



3(x- 1) 



In 2 x 



(17) 



Gf(x) 



^(3), 



+ A« 



A 



(2) 



21 _4«(8 + 4«) 



61-39x (7-60x-14x 2 



12 (x- If 



6 (x 



In x — 



3(x- 1) 
14 x 



3 (x 



In x 



X 



1?; 



^(F 8 (1) (x)-3F 7 (1) ( 



31n^ + Aj? - 1 



32 
27 



x , 



Af(x) = ^(Ff (x)-3Ff(x; 



x 



28 pW, 



X 



x 



3 8 



x 



(19) 

(20) 
(21) 



The various functions A appearing in eqs. (|15|)-(|19|) contain the effect of the renor- 
malization of the chargino-stop-quark vertex due to the 0(fi g ) heavy particles. Their 
explicit expressions are given in appendix A. In the same equations the terms not 

(1 3) 

proportional to the one-loop functions F 7 g represent the /i-independent part of the 
contributions coming from two-loop squark-chargino diagrams. The [i w dependence 
in eq. (|13|) and eq. ([14]) satisfies the relation 



1 m - dC\ 0) 



+ 

z i=i 



Ml 



(P)eff c (0) 



2 

£ 



m 



1 2 



16 



V 2 cos p m Xj K 



(22) 



that ensures that physical observables are independent of \i w to 0(a s ). In eq. ( p2| ) 
7™ = 8 is the LO anomalous dimension of the top mass. 

Finally, we report the effects of O(fig) heavy particles in diagrams involving W, 
physical and unphysical charged scalars exchanges. As previously discussed, in our 
approximation they are introduced as a renormalization of the relevant couplings 
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and the corresponding contributions are given by 

^Cf e// GM = ^ (W t b + W t s ) G? (tw) - P (W b c + W°) (23) 
S w cP eff (f, w ) = ^(w' + W t s )G^(tw)-^(w b c +W:) (24) 
8*Ci¥"M = {HI + H*) F$ (th) 

+ ^(H s t + H b ) F§ (th) (25) 
hMVefff.. \ - 4 /'rr* i z?W C+„.A _ 4 m* i f( 2 ) 



2 Q; 8 ;e// (M = g^ + ^j^M-g^ + ^nsM • (26) 
Here tw = m t (fi w ) 2 /M^, th = m t (^ w ) 2 /M|±, 

. 23 - 67a; + 50x 2 x(2 - 7x + 6x 2 ) , 
G 7 (x) = — 1 — ; lnx; 



' T 1 ' 36(x-l) 3 6(x- l) 4 

- 5x — 5x 2 x(l — 2q 
12(x-l) 3 + 2(x-l) 



„m, , 4 — 5x — 5x 2 x(l — 2x), , . 

Gf* = 7777 777- + 777 77J ln ^ ( 27 ) 



the functions F^g are defined in eqs. (54)-(55) of ref. |7j], and the renormalization 
functions W, H, U are given in appendix A. 



4 Numerical Analysis 

In this section we present a numerical analysis of BK(B — > X s j), with special 
attention to the comparison between LO and NLO results in the supersymmetric 
model. We focus on the situation in which our approximation is adequate, i.e. light 
stop and charginos, and show how the strong bounds on the charged Higgs mass 
can be relaxed. 

We use the same numerical inputs and conventions as in ref. 0, except for 
the photon-energy cut off F 7 > (1 — 5)m&/2, which is now chosen to be 5 = 0.90 
(see discussion in ref. [[HJ). In our numerical analysis we also include the leading 
logarithmic electromagnetic effects as computed in refs. |], |TU|. They amount to a 
decrease of the calculated branching ratio of about 7% to 8%. We neglect all other 
two-loop electroweak effects which appear to be very small p? . 



Light supersymmetric particles can affect the measured values of the B° — B° 
mixing AMg and of the CP-violating ex parameter and, ultimately, the extraction 
of the CKM angles, which enter the calculation of BK(B — > X/y). Within our 
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assumptions, it is not difficult to compute this effect. The dominant supersymmetric 
contributions to AMb and ex come from box diagrams with chargino/light stop 
and charged Higgs/top quark exchange. Under the assumption of minimal flavour 
violation, these contributions have the same structure of CKM angles as the top- 
quark box diagram in the SM. Since the dominant effects on AMb and ex in the 
SM indeed come from the top quark, the CKM-angle dependence is the same for all 
contributions, and the effect of supersymmetry is parametrized by 



R = l + 



A 



SM 



(2f 



Here Asm is the amplitude of the top-quark box diagram, and Asusy is the ampli- 
tude of the chargino/light stop and charged Higgs/top quark box diagrams, whose 
analytical expressions are given in the appendix B. The crucial point is that R is 
independent of CKM angles. In the whole supersymmetric parameter space, R turns 
out to be larger than one. 

In terms of the Wolfenstein parameters rj and p, the CKM structure of AMb 
and 6k in supersymmetry is 

~:i-p) 2 

77(1- 



VuV tb \ 2 R 



e K ~ \hn{V*VuY\ R 



+ 7] 

p)R. 



R 



(29) 
(30) 



The new physics contribution can be effectively reabsorbed in the Wolfenstein pa- 
rameters (1 — p)y/~R and rjvR. This means that the "true" values of rj and p are 
related to the values rj SM and psm extracted from the usual SM fit by 

VSM 



V 



p=l 



R 

(1 - Psm) 



R 



(31) 



(32) 



We have checked that eqs. (|3T|) — d32|) give an extremely good approximation of a 
complete numerical fit. The combination of the CKM angles entering the evaluation 
of BR(B -> X s7 ) is 



X = 



v t *v tb 



cb 



1 + (2p - 1)A 2 = 1 



1 



'1 - Psm) 



(33) 



where A = \V US \ does not depend on the supersymmetric masses. Therefore, in the 
presence of supersymmetry, the value of X is related to the usual SM input X SM by 

1 \ 



X = X SM + (A + 1 - X 



SM 



1 



R, 



(34) 
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This is the value of X to be used in a supersymmetric analysis of BK(B — > X s ~f). 

It is known that the chargino/stop contribution to the Wilson coefficient Cj can 
interfere with the SM and charged Higgs contributions either constructively or de- 
structively, depending on the parameters choice. The case of destructive interference 
is particularly interesting, since it is possible to make a light charged Higgs consistent 
with the measurement of BR(I? — > X s r {) . For illustration, we choose the supersym- 
metric parameters corresponding to a strong cancellation between the chargino and 
charged-Higgs contributions, and plot in fig. 1 the predicted BR(B — > Xgj) as a 
function of fi g . The LO result in supersymmetry is compared with the result ob- 
tained in the limit \i g — > oo (denoted by LOdecoupled in fig. 1). This shows how the 
asymptotic behaviour is reached, as \x g grows. The NLO result is also shown in fig. 1. 
In this case, we cannot define an asymptotic behaviour, as the NLO branching ratio 
hardly decreases with \i g and actually starts rising again for ji g > 0.7 TeV. This 
non-decoupling effect is not surprising pB| . It appears because, at the NLO, the 
chargino vertex contains terms proportional to log(/ig//i w ). Indeed, the chargino 
coupling constant is related to the ordinary gauge and Yukawa couplings only by 
supersymmetry. In the effective theory below fi g , supersymmetry is explicitly bro- 
ken and the renormalization of the chargino vertex develops a logarithm of the large 
supersymmetry-breaking scale. Notice the important effect of the NLO corrections 
calculated here. For the supersymmetric parameters of fig. 1 and for /i§ = 1 TeV, 
the LO branching ratio is 1.9 x 10~ 4 , while the NLO result is 3.7 x 10~ 4 . Had we ne- 
glected in the NLO analysis our new supersymmetric contributions to the matching 
conditions, the result would have been 2.4 x 10~ 4 . The large shift induced by the 
NLO corrections to Cjj/(fi w ) is mainly due by the large non-decoupling logarithms 
of fj, g /fi w . Finally, in fig. 1 we also show the LO and NLO results obtained in a two- 
Higgs doublet model with charged-Higgs mass equal to the supersymmetric case. 
The light charged Higgs (taken with a mass of 100 GeV in fig. 1) is clearly incom- 
patible with the measurement of BR(B — » X s j), in the absence of an appropriate 
chargino contribution. 

In order to assess the impact of the different contributions, it is also interesting 
to consider the NLO Wilson coefficients C^J at the weak scale \i w = Mw in the 
supersymmetric configuration chosen in fig. 1 with all heavy masses set to 1 TeV, 
and compare them with the corresponding quantities in the SM and the two-Higgs 
doublet model. The results are shown in the table 1. Notice that, in the super- 
symmetric case, the QCD corrections to Cy amount to about 60%. Such a large 
correction is clearly related to the fact that the LO approximate cancellation among 
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different contributions is partially spoiled at the NLO level. Indeed, in the case un- 
der consideration, the NLO effects increase the SM contribution to \C^ f {M w )\ by 
almost 10% and decrease the charged-Higgs and the chargino contributions by about 
20% and 30%, respectively. This situation of partial cancellation and enhanced sen- 
sitivity to NLO correction is actually the most interesting phenomenologically. It 
is under this condition that a light charged-Higgs mass is still allowed and that our 
calculation is essential. 





C e 7 ff (M w ) 


C e 8 ff (M w ) 


LO SM 


-0.198 


-0.098 


NLO SM 


-0.220 


-0.119 


LO 2HDM 


-0.529 


-0.336 


NLO 2HDM 


-0.493 


-0.326 


LO SUSY 


-0.143 


-0.141 


NLO SUSY 


-0.229 


-0.183 



Table 1: Wilson coefficients evaluated at the scale M w in the Standard Model (SM), 
the two-Higgs doublet model (2HDM) and supersymmetry (SUSY) at LO and at NLO 
for tan/3 = 1, M H ± = m h = m X2 = 100 GeV, m xi = 300 GeV, Q~ t = -vr/10, A b = A t ; 
all other squarks and the gluino are degenerate with mass of 1 TeV. The lighter chargino 
is predominantly higgsino. 

We now want to investigate the region of parameter space in which we can 
significantly relax the bound on the charged Higgs mass from BR(£> — > X s ^). We 
have scanned over the relevant supersymmetric parameters, assuming a stop mixing 
angle less than 7r/10 in absolute value, consistently with our assumption of a mainly 
right-handed light stop. In fig. 2 we show, for tan (3 = 2 and 4, the maximum values 
of the lighter chargino and lighter stop masses for which a charged Higgs of 100 
GeV is consistent with the 95% C.L. CLEO limit of 4.2 x 10" 4 . We have chosen 
this value to give a conservative estimate. With improved experimental results and 
revised analyses of the photon energy spectrum, this bound may become significantly 
more restrictive. The effect of the NLO corrections is to make the upper bounds 
on the chargino and stop masses quite more stringent. In the same figure we also 
present the results obtained using the renormalization group evolution and the SM 
and charged-Higgs contributions to the Wilson coefficients at the NLO, while the 
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purely supersymmetric effects on C?f 8 (/jl w ) are evaluated at the LO. This is to show 
the impact of the calculation presented here. As seen from the figure, the effect 
of the QCD corrections to the supersymmetric contribution to Off/ (f/, w ) is very 
significant. 

5 Conclusions 

In this paper we have computed the QCD NLO corrections to the matching con- 
ditions of the Wilson coefficients relevant for BR(B — ► X s ^). We have used some 
theoretical assumptions to simplify the result and to concentrate on the most rele- 
vant part of the supersymmetric parameter space. We have assumed minimal flavour 
violation or, in other words, that the flavour-violating interactions in supersymme- 
try are primarily in the charged- current sector and are determined by the CKM 
angles. This is often a good approximation in a variety of models. At worst, our 
result represents an unavoidable contribution to be added to new effects coming 
from other sources of flavour violation. This assumption allow us also to limit the 
number of unknown parameters. 

We have focused on the case in which the purely supersymmetric contribution 
is sizeable and can compensate the effect of a light charged-Higgs boson. For this 
reason, we have assumed that charginos and a mainly right-handed stop are consid- 
erably lighter than the other squarks and the gluino. This assumption has allowed us 
to use an effective theory, in which the heavy particles are integrated out, retaining 
only the first term in an expansion in the ratio between light and heavy masses. 

In sect. 3 we have presented the analytic results for the NLO matching conditions 
of the Wilson coefficients in supersymmetry, under the approximations stated above. 
We have also performed a numerical investigation of the supersymmetric parameter 
region in which there is a sizeable cancellation among various different terms. In 
this region the impact of the NLO corrections to the supersymmetric contribution to 
the Wilson coefficients is very significant because of two conspiring effects. On one 
side there is a large renormalization of the one-loop supersymmetric contribution, 
mostly coming from logarithms of the ratio between the high and the intermediate 
mass scales. On the other side, the effect of the large NLO corrections to the 
supersymmetric contribution is greatly enhanced whenever there is an approximate 
cancellation at the LO. 
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Appendix A 

In this appendix we present the results for the renormalized vertices in the effective 
theory at scales between \i~ g and \x w . In deriving the various corrections we exploited 
the fact that, in our framework, the light stop is mainly right-handed and therefore 
we can neglect terms O (sin 6 /i w / )j,g) . However in the formulae below we do not 
expand functions of m-^/m- and we keep the explicit sin 2 Q~ t terms, because we feel 
it will be easier for the reader to understand the origin of the different terms. 



w 

The renormalized du W vertex is 
where 

W d = 

2 

with x\ = nr-^/rri-, Wi = rn\./m~ and 
W[x,y] 



= \ (cos 2 9u W[ Xl , w x \ + sin 2 ^ W[x 1: 



w 2 



v + y-2xy x 3 - 2xy + x 2 y 2xy - xy 2 - y 3 

In x + 77777 7 In y. 



(x-l)(y-l) (x-l) 2 (x-y) ( y _ 1)2(3 _ j,) 

For equal masses W[x, x] = 0. 



Unphysical scalar 

The renormalized dtd> + vertex is 



(-0 



9 



td 



where 



U d = U H^x,} - cos 2 OiH^u,}- sin 2 OiH^ 



n, <: 1 • y^ r '') " h "- ( ' • 'i^ 1 ' 
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-2 Ad -^^ H 2 [ Xl , X2 ] 



m g 



+2 Ad * tan/3 (cos 2 g t -g 2 [ Ul , x 2 ] + sin 2 6~ t H 2 [u 2 , x 2 }) , 



m-g 



Uf = -^(^[xil-cos^f^tuxl-sin 2 ^-^^]) 



2 

A, - // cot /? 

-2 if 2 K> U 2 

A u — cot /3 



+2 ^ ^""^ (cos 2 t ~tf 2 K xi] + sin 2 % tf 2 K 
m 5 



with x 2 = m~^/m~, U; L = m~ /m~ and 



rr r 1 1 2X-X 2 , 

#iF = i + 71 « ln x 

1 — x (1 — re)" 1 



#2F,y] = ^ rlnx + - ^ rlny 

(l-x)(x-y) {l-y){y-x) 

Physical scalar 

We write the renormalized vertex d t h + as 

72% ytd [ tan 0mda -{ 1 + ^ Hd ) + cot ^ m < a+ i 1 + £ ^ 



where 



H d = ^(H 1 [x 1 ]-cos 2 e i H 1 [u 1 ]-sm 2 e i H 1 [u 2 \) 

A -//tan/3 
! H 2 [x 1 ,x 2 \ 

nm - 



-2 



""3 

+2 A, + p cot g , 2 u + gin2 v 

^ = —(^[Xil-cos^f^tuxl-siii 2 ^^^]) 

o Au-H cot/3 

-2 # 2 [ui,u 2 \ 

m- g 

+2 A u + f i tan/3 , 2 + s , n2 u v 
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Chargino 

Unlike the case of the W and of physical and unphysical scalars, in the chargino 
sector the 0(a s ) corrections due to gluon and gluino exchanges are not separately 
finite. Therefore, as shown below, the renormalization of the chargino vertex due 
to the 0(ng) heavy particles will contain 1/e poles (e = (4 — n)/2, with n the 
dimension of the space-time). The cancellation of these poles against similar terms 
coming from two-loop diagrams in which a gluon is present provides a check of our 
calculation. We write the corrected Xa d t 2 vertex as 

(—iq) V t d \ — sin Qz U a2 — f= a+ ( 1 H C c 



+ 

with (ft is the 't-Hooft mass 



V2M w smf3 a - V + 3vr ' ' V ' 3vr 



cos 9- t V a2 . a h+-J-CZ d )+ sin 9- t V al a_ 1 + -a- CL 



(Al) 



3 3, (m 2 h 



c " - 4 + h#) +as+% - (A2) 

The functions in eq. (|A2j), also appearing in eqs. (p~5|) — (P~9|) , are given by 

A?J = - 

A ?' = ^+(l-oot%^i)/f,[x 2 ] + iH 1 [x a ] + iH ] [r 2 | 

_ 2 A d - M tan^ 
'"v 

» f2i 5 r 1 Tr r , 1 TT r , A u — u, cot 3 TT r , 
A! J = - + // 3 ki] + -^iN + 2 i/ iN-2^^ -iT 2 [«i,w a ] 




'a 



where 



H 3 [x] = — — In x. 
1 — x 



In eq. ( |A2|) the factors r\ Y = —3/2 and r] g = —1/2 are induced by the fact that the 
MS renormalization does not preserve supersymmetry, as discussed in the text. 
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Appendix B 



This appendix contains the expression of the parameter R defined in eq. (0), which 
can be derived from the results of ref . |16[ . The SM contribution is 

tw 3 - Yltw 2 + 15 tw - 4 + Qtw 2 \ntw .„ . 

A '» = U^if < B1 ) 

with tw = m t (fi w ) 2 /M^. The supersymmetric contribution can be split into charged 
Higgs and chargino contributions: 



^SUSY 



A H + A. (B2) 



a Ath 2 - 1 - 2 th \nth)th , „ 

A H = cot 4 ft — — —: — — + 2 cot 2 tw 

4 [tn — I) 6 



F'(tw,hw) + -G'(tw,hw) 



(B3) 



with th = m t (n w ) 2 /M 2 ± and hw = M 2 ± /M^. The functions F' and G' are given 
by 



F'{x,y) 
G\x,y) 



y in y 



(x 2 — y) lnx 

(x - y) 2 (x - l) 2 (x-y) 2 (y-l) (x-y)(x-l) 



x 



(x - y){x - 1) 



i y 
+ 



x — 1 x — y 



\nx 



y 2 In j/ 



(x-y) 2 (y- 1) 



(B4) 



Consistently with our assumptions, we keep only the light stop contributions to A: 



M 4 

^2^2 C 2i r 2j ^ U/2j) ^ij;, 
,3=1,2 m t m Xj 



(B5) 



where 7/ 2 j is the variable defined in eq. (|7|) and is the ratio of the squared masses 
of two charginos. 
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Figure 1: Dependence of the prediction for the B — > X s ^ branching ratio on the 
heavy mass scale fig for the following choice of parameters: tan/3 = 1, M H ± = mi 2 = 
m X2 = 100 GeV, m Xl = 300 GeV, 6* t ~ = — 7r/10, Ab = A t , all heavy particle masses 
equal to fi g ; the lighter chargino is predominantly higgsino. The value corresponding to 
" LOdecoupled" is the result of the LO calculation in the limit in which the heavy masses 
decouple. The two upper lines are the LO and NLO predictions in the two-Higgs doublet 
model (2HDM). 
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Figure 2: Upper bounds on the lighter chargino and stop masses from the CLEO 95% 
CL limit on BR(B -> X s -f) in the case M H ± = 100 GeV. We have taken \6 t ~\ < tt/10, 
\fji\ < 500 GeV, At, = A t , and set all heavy masses to 1 TeV. For tan/3 = 2 and 4 we 
show the results of the LO and NLO calculations. The result of neglecting the new NLO 
supersymmetric contributions to the Wilson coefficients is labelled as "NLO running". 
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